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In this report we have analyzed a simple effective model for a description of magnetically ordered 
insulators. The Hamiltonian considered consists of the effective on-site interaction {U) and the 
intersite Ising-like magnetic exchange interaction (J) between nearest neighbors. For the first time 
the phase diagrams of this model have been determined within Monte Carlo simulation on 2D-square 
lattice. They have boon compared with results obtained within variational approach, which treats 
the on-site term exactly and the intersite interactions within mcan-fiold approximation. We show 
within both approaches that, depending on the values of interaction parameters and the electron 
concentration, the system can exhibit not only homogeneous phases: (anti-)ferromagnetic (F) and 
nonordered (NO), but also phase separated states (PS: F-NO). 

PACS numbers: 71.10.Fd, 75.30.Fv, 64.75.Gh, 71.10.Hf 



I. INTRODUCTION 

The extended Hubbard model with spin exchange in- 
teraction [1-4] is a conceptually simple effective model 
for a description of magnetically ordered insulators in 
narrow band systems. 

In this report we will focus on the zero-bandwidth limit 
of the extended Hubbard model with magnetic Ising-like 
interactions, which has the following form: 

H = U^ni^hii^-2J^sm- iJL^hi, (1) 

where U is the on-site density interaction, J is z- 
component of the intersite magncitic exchange interac- 
tion, XI (ij) restricts the summation to nearest neigh- 
bors, denotes the creation operator of an electron 
with spin a at the site i, hi = iiia, fii„ = cf^Cia and 
— K'^it ~ ^ii)- The chemical potential fj, depending 
on the concentration of electrons is calculated from 

i 

with < n < 2 and N is the total number of lattice sites. 

The hamiltonian (1) can be considered as a very simpli- 
fied model for the family of A0.5M2X4 compounds (where 
A is Ga or Al, M is one of the transition metals V or 
Mo, and X is S, Se, or Te). These compounds exhibit 
very interesting ferromagnetic behavior which is a mix- 
ture of itinerant and localized behavior [5, 6]. However, 
the single-particle excitations play a dominant role in the 
magnetic behavior of these compounds. Although the 
electrons arc not itinerant in the system, there is a finite 
density of states at the Fermi level, and therefore low en- 
ergy charge excitations are possible [7] . This points that 
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the magnetic properties result from this band of localized 

electrons. 

The model (1) can be treated as an effective model of 
magnetically ordered insulators. The interactions U and 
J will be treated as effective ones and be assumed to in- 
clude all the possible contributions and renormalizations 
like those coming from the strong electron-phonon cou- 
pling or from the coupling between electrons and other 
electronic subsystems in solid or chemical complexes. In 
such a general case arbitrary values and signs of U arc 
important to consider. One should notice that ferromag- 
netic ( J > 0) interactions are simply mapped onto the 
antiferromagnetic cases (J < 0) by redefining the spin di- 
rection on one sublattice in lattices decomposed into two 
interpenetrating sublattices. Thus, we restrict ourselves 
to the case J > 0. 

For the model (1) only the ground state phase dia- 
gram as a function of /i [8] and special cases of half-filling 
(n = 1) [9] have been investigated till now. Some our pre- 
liminary results have been also presented in Ref. 10. 

Wc have performed extensive study of the phase dia- 
grams of the model (1) for arbitrary n and /U [10-12]. In 
the analysis we have adopted two complementary meth- 
ods: (i) a variational approach (VA), which treats the on- 
site interaction term (U) exactly and the intersite inter- 
actions (J) within the mean- field approximation (MFA) 
and (ii) Monte Carlo (MC) simulations for d = 2 dimen- 
sional square (SQ) lattice in the grand canonical ensem- 
ble. In this report we present some results concerning 
strong on-site repulsion limit. 

The ferromagnetic (F) phase is characterized by non- 
zero value of the magnetic order parameter (magnetiza- 
tion) defined as m = (1/-^) X]j(sf)- In the nonordered 
(NO) phase m = 0. 

The phase separation (PS) is a state in which two do- 
mains with different electron concentration exist (coexis- 
tence of two homogeneous phases). In the model consid- 
ered only one PS state (F-NO) can occur. 

In the paper we have used the following conven- 
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FIG. 1. Phase diagrams ksT j vs. /i/Jo for ?7/Jo = 1, 10 
(as labeled) obtained within VA. Dotted and solid lines in- 
dicate first and second order boundaries, respectively. T de- 
notes a tricritical point. 

tion. A second (first) order transition is a transition 
between homogeneous phases with a (dis-)continuous 
change of the order parameter at the transition temper- 
ature. A transition between homogeneous phase and PS 
state is symbolically named as a "third order" transition. 
During this transition a size of one domain in the PS state 
decreases continuously to zero at the transition temper- 
ature. One should notice that the first order transition 
line on the diagrams for fixed [i splits into two "third or- 
der" lines and it is connected with occurrence of PS sates 
on the diagrams for fixed n. 

We also introduce the following denotation: Jq = ^ 
where z is the number of nearest neighbors. 

Obtained phase diagrams are symmetric with respect 
to half-filling (n = 1, /i = /i — C//2 = 0) because of the 
particle- hole symmetry of the Hamiltonian (1), so the 
diagrams will be presented only in the range < rt < 1. 



II. RESULTS AND DISCUSSION 
A. The variational approach 

In this subsection we discuss the results for strong 
on-site repulsion obtained within VA. The dependencies 
of the transition temperature F-NO as a function of /l 
for [// Jo = 1 and C// Jo = 10 (this is close to the limit 
[/ — )■ -l-oo) are shown in Fig.l. The range of the F phase 
stability is reduced with decreasing [//Jo- A tricritical 
point T connected with a change of the F-NO transition 
order is located at fcsT/Jo = 1/3 and its location does 
not dependent on U j Jq in the limit considered. 

If the system is analyzed for fixed n [10], at sufficiently 
low temperatures the homogeneous phases are not states 
with the lowest free energy and the PS state can oc- 
cur. On the phase diagrams, there is a second order 
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FIG. 2. Phase diagrams fcsT/Jo vs. n for: (a) U/Jo — 1 and 
(b) U/Jo = 10 obtained within VA. Solid and dashed lines 
indicate second order and "third order" boundaries, respec- 
tively. T denotes a tricritical point. 



line at high temperatures, separating F and NO phases. 
A "third order" transition takes place at lower tempera- 
tures, leading to a PS of the F and NO phases. The crit- 
ical point for the phase separation (T) lies on the second 
order fine F-NO and it is located at ksT/Jo = 1/3 and 
n= 1/3. Phase diagrams for U/Jq = 1 and U/Jq = 10 
are shown in Fig. 2. With increasing fc^T/ Jq the system 
exhibits either a sequence of transitions: PS— >-F— ^NO 
(for 1/3 < n < 1) or a single transition: PS— >N0 (for 
n < 1/3) and F^NO (for n = 1). 



B. Monte Carlo results 

Here, we present a numerical investigations of model 
(1), using standard MC methods in the grand canonical 
ensemble (for details see e.g. Ref. 13). The MC sim- 
ulations have been done for two dimensional SQ lattice 
(z ~ 4) with periodic boundary conditions. The size of 
the lattice is relatively small, i.e. 10 x 10. 

The general properties of MC phase diagrams are sim- 
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FIG. 3. Phase diagrams ksT / Jo vs. Jl/Jo for U / Jo = 1, 10 
(as labeled) obtained within MC simulation for 10 x 10 SQ 
lattice. T denotes a tricritical point. Details in text. 



ilar to those obtained within VA. However, it is obvi- 
ous that the transition temperatures resulting from MC 
simulations for d = 2 SQ lattice are lower than those ob- 
tained within VA, which is exact in the limit of infinite 
dimensions. 

The phase diagrams as a function of /l for C//Jo = 1 
and U / = 10 are shown in Fig. 3. The transitions in fi- 
nite systems are not sharp (the finite-size effect on the or- 
der parameter, i.e. in the NO phase m ^ is larger than 
zero near the transition) and thus the precise location of 
boundaries between different phases are determined by 
the discontinuity of the magnetic susceptibility. A tri- 
critical point T connected with a change of the F-NO 
transition order is located at ksT/Jo sa 0.27. The F-NO 
transition can be first order (for temperatures below 
T-point) as well as second order (for temperatures above 
T-point). The maximum for the F-NO transition tem- 
perature is at half- filling (n ~ 1, fi/ Jq = 0) and it equals 
(i) fcsT/Jo « 0.48 for U/Jq = 1 and (ii) fcsT/Jo « 0.59 
for U/Ja = 10. The behaviors of the boundaries at low 
temperatures (i.e. for kgT/jQ < 0.1) have not been de- 
termined. 

One can translated the (grand-canonical) diagrams 
from Fig. 3 into the (canonical) diagrams for arbitrary n 
by the standard way. The resulting diagrams are shown 
in Fig. 4. At higher temperatures the F and NO phases 
are separated by a second order line. At lower tem- 
peratures (below T-point) the PS state occurs, which 
is separated from homogeneous phases (i.e. F and NO 
phases) by "third order" boundaries. The tricritical point 
is placed at substantially higher electron concentrations 
(n « 0.61) in comparison to VA results and (as in VA) 
its location is independent of the on-site repulsion U/Jq 
in the limit considered. 
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FIG. 4. Phase diagrams kBT/Jo vs. n for: (a) U/Jo — 1 and 
(b) U/Jo = 10 resulting from MC simulation for 10 x 10 SQ 
lattice. Denotations as in Fig. 2. Details in text. 



III. FINAL COMMENTS 

We considered a simple model of magnetically ordered 
insulators. We presented phase diagrams for strong on- 
site repulsion including a tricritical behavior obtained by 
Monte Carlo simulations and compared them with VA 
results. It was shown that MC results are qualitatively 
similar to those derived within the VA. However, one 
should notice that the MC transition temperatures are 
significantly smaller than VA ones. The F-NO transi- 
tion can be second as well as first order. At sufficiently 
low temperatures, where the F-NO transition is discon- 
tinuous (if P/Jq is fixed), homogeneous phases do not 
exist (if n is fixed) and the phase separated states have 
a lowest energy. 

Let us stress that the knowledge of the zero-bandwidth 
limit can be used as starting point for a perturbation 
expansion in powers of the hopping and as an important 
test for various approximate approaches analyzing the 
corresponding finite bandwidth models. 

We leave the problem of detailed analysis concerning 
arbitrary U/Jq for the future investigations [12]. 
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